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ativistic current 
2] for electronic 



Among the extensions of density functional theory (DFT), the non-re 
density functional theory (CDFT), designed by Vignale and Rasolt (VR) 
systems in a static magnetic field, stands out for its conceptual subtlety. Electrons couple to 
a magnetic vector potential A(r) via the so-called paramagnetic current density j p (r), which 
in the non-relativistic theory differs from the physical current density j(r) and is related to 
the latter by 

j(r) =j p (r) + p(r)A(r) (1) 



[Here, „(r) is the gronnd-state density, and we nse units in which e 



m 



1] . It has been 



proved [l| that p(r) and j p (r) uniquely determine (through a one-to-one correspondence) the 
ground-state wave function, but not necessarily the external scalar and vector potentials, 
V(r) and A(r), that give rise to them. Schematically, one has 



(2) 

potentials are not uniquely 



VR: {p,j p }^^^{V,A}, 
where the last arrow points only one way. Although the externa 

determined by the densities - a fact that was overlooked in Ref . [l and was first pointed out by 
Capelle and Vignale in Ref. O - the one-to-one correspondence between densities and ground- 
state wave functions is sufficient to construct the universal "internal energy" functional 
F[p, j p ] and the exchange-correlation energy functional E xc p, j p ], and the constrained search 



algorithm and the variational principle remain unaffected 



4] 



In spite of this, a recent paper by 



'an and Sahni (PS) 



5| seizes on the non-uniqueness 



of the potentials - the point of Ref. [3| - as a casus belli against the VR version of CDFT. 
The authors of this paper have stated that (i) the ground-state wave function "cannot be 
a functional of p and j p " (end of Section 6) and (ii) the physical current j is the "natural 



basic variable" for CDFT 



6]. 



Statement (i) is obviously a mistake: the ground-state wave function is indeed a unique 
functional of p and j p [see Eq. (J2J) above] , as demonstrated in Refs. In this comment 

we show that statement (ii) is also erroneous. The physical-current density version of CDFT 
leads to problems far more serious than the ones it purports to cure. Let us see why. 

The essential difficulty with the choice of the physical current j as a basic variable is 
that the correspondence between density and wave functions, which is one-to-one in the VR 
formulation, becomes many-to-one in the PS formulation. Schematically, PS have proved 



PS: { p ,j}^{y )A }^ 



(3) 



where the last arrow points only one way. We believe this statement to be correct, since it 



follows immediately from previously proved results js]. However, different current densities 
may now be associated with one and the same ground-state wave function. Indeed, given 
a wave function, it is impossible to work out the physical current without the additional 
knowledge of the vector potential. The consequences of this fact are disturbing, as we shall 
now show. 

Let us follow PS and write the energy functional as 

EvAp'A'] = Fp S [p',y] + J j'(r) ■ A(r)dr -I J p\v)A 2 {v)dv + j p'{v)V{v)dv , (4) 

where Fp$[p,j] is the expectation value of the standard many-body Hamiltonian T + U = 
~\ J2i ^l~^~J2i<j U{vi — Tj) in the ground-state wave function, which is uniquely determined 
by n(r) and j(r) according to ([3]). We note that, contrary to what is stated in Eq. (59) 
of PS, the functional Fps[p',j'} cannot be constructed from a constrained minimization of 
T + U on the set of wave functions that yield p and j. The reason is that this set is undefined, 



unless the vector potential is specified - and therefore Eq. (59) of Ref. 5I is meaningless. 

The standard variational principle for ground-state wave functions mandates that we 
determine the ground-state values of p and j by minimizing Ey,A[p',y] with respect to p' 
and j' at constant V and A. Notice, however, that j = j p + pA [Eq. (JTJ)]. From the form 
of this relation we clearly see that a minimization with respect to j and p at constant A 
and V is equivalent to a minimization with respect to j p and p under the same conditions. 
Once A is fixed, j and j p become essential the same variable, and we may as well rewrite 
the functional to be minimized in terms of j p as follows: 

E VA [P', Jp] = F VR [p' : y p ] + J j;(r) • A(r)dr + I J p\r)A 2 {v)dv + J p'{v)V{v)dv (5) 

(notice the change in sign of the term containing A 2 ). Fyji[p',y ] is precisely the universal 
VR energy functional of j p . Thus, the minimization of the PS functional with respect to j 
at constant A is equivalent to the minimization of the VR functional with respect to j p at 
constant A, if one takes proper care of the A that is hidden in the definition of j. Under 
the stipulation of constant A, the PS theory is nothing more than a garbled reformulation 
of the VR theory. 

However, in the above discussion we took a lenient view of the functional (Jl|). A stricter 
interpretation, as advocated by PS, would be to treat j as the "natural basic variable" and 



perform the minimization with respect to it, only keeping constant the A that explicitly 
appears in Eq. fll]). By doing this, we run into a serious problem. We know from Ref. [3 
that it is possible to find two different sets of infinitesimally close densities, {p(r), j(r)} and 
{p(r), j'(r) = j(r) + <5j(r)}, that correspond to the same ground-state wave function, and, 
therefore, yield the same F: Fps[p,j] = Fps[p,j']. Then in the infinitesimal variation from 
{p(r), j(r)} to {p(r), j'(r)} at constant V and A we have 5Fp$ = and bp = and the full 
energy functional changes by 

6E ViA = J 5j(r)-A(r)dr^0. (6) 

The energy functional is no longer stationary for infinitesimal variations about the chosen 
ground-state densities p and j! This completely spurious result has its root in the ignorance 
of the vector potential-dependence of j and indeed disappears as soon as one rewrites the 
energy functional in the equivalent VR form of Eq. (jSJ), for which SFvr and 5j p are both 
zero. But, the paradox shows how dangerous it is to work with a basic variable that is not 
an intrinsic property of the system, and therefore cannot be uniquely determined by the 
wave function alone. 

To summarize, we have found that the PS formulation of CDFT is, at best, a garbled 
reformulation of the VR theory, and, at worst, a potential source of mistakes insofar as it 
complicates the formulation of the variational principle and prevents the constrained search 
construction of the universal functional. 
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have j p = j' and p = p' , and since the vector potentials must necessarily be different [3[ we 
immediately conclude that j 7^ j'. 
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